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Abstract. The rapid rotation of the planets causes cyclonic 
thermal turbulence in the cores, which can be responsible for 
the generation of the large-scale magnetic fields observed out 
of the planets. We consider the model which lets us to repro- 
duce the typical features of the small-scale geostrophic flows 
in the physical and wave spaces. We present estimates of ki- 
netic and magnetic energies fluxes as a function of the wave 
number. The joined existence of the forward and inverse cas- 
cades are demonstrated. We also consider mechanism of the 
magnetic field saturation at the end of the kinematic dynamo 
regime. 



1 Introduction 

Many astrophysical objects such as galaxies, stars, the earth, 
and some of the planets have large scale magnetic fields 
that are believed to be generated by a common universal 
mechanism - the conversion of kinetic energy into mag- 
netic energy in a turbulent rotating shell. The details, how- 
ever, - and thus the nature of the resulting field - dif- 
fer g reatly. The challenge for th e dynamo theory see, 
e.g., JHoUerbach and Riidigen, 120041) . is to provide a model 
that can explain the visible features of the field with re- 
alistic assumptions on the model parameters. Calcula- 
tions fo r the entire plane t are done either with spectral 
models ( Kono and Roberts, 2002) o r finite volume methods 
(iHejda and Reshetnya k. 2004; Harden and HansenL l2005b 
and finite difference CKagevama and Sato|, 1 1997b . and have 
demonstrated beyond reasonable doubt that the turbulent 3D 
convection of the conductive fluid can generate a large scale 
magnetic field similar to the one observed out of small ran- 
dom fluctuations. However, both these methods cannot cover 
the enormous span of scales required for a realistic param- 



eter set. Even for the geodynamo (which is quite a mod- 
est on astrophysical scales case) the hydrodynamic Reynolds 
number estimated on the west drift velocity is Re ^ 10^. 
In addition planets are the rapidly rotating bodies. So, the 
time scale of the large-scale convection in the Earth's core 
is ^ 10"^ years, during which the planet itself makes ^ 10^ 
revolutions (in other words the Rossby number Ro ~ 10^^). 
As a result, there is an additional spatial sc ale ^ E"'^^'^ L, 
where E ~ lO"'^^ is the Ekman number (IChandrasekhan 
1961; Bu sselll970h . associated with the cyclonic structures 
elongated along axis of rotation, which is quite larger then 
the Kolmogorov's dissipation scale Id ^ Rc^ ' L however 
which is still too small to be resolved in the numerical simu- 
lations with present resolution / ^ (10^^ ~ lG^^)i with L 
for the large scale. 

Presence of the rapid rotation leads not only to the 
change of the spatial uniform isotropic Kolmogorov like 
solution to the quasi-geostrophic (magnetostrophic) form 
but to rather more fundamental consequences. Rapid rota- 
tion leads to degeneration of the third dimension (along the 
axis of rotation) and can cause inverse cascade in the sys- 
tem. The inverse cascades is a well-known phenomenon in 
the two-dimensional turbulence and is a good example of 
self-organization, when the large-scale structures are feeded 
by the small-scales tur bulence (Kraichnan and Montgomery! 
1980t iTabelingl 120021) . So far, that the quasi-geostrophic 
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state formally still is a three dimensional (the vector fields 
have three components) behavior of such systems can differ 
from the two-dimensional flows as well as from the three di- 
mensional ones. It can happen, that quasi-geostrophic turbu- 
lence can exhibit simultaneously features similar to the both 
extreme cases: 2D and 3D. Further we consider behavior of 
the fluxes of the energy in the wave space for the well-known 
in geodynamo regimes based on the Boussinesque thermal 
convection. For simplicity we consider the Cartesian geom- 
etry which is simpler for modeling of the rapidly rotating 
dynamo systems and was used in many researches on geo- 
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dynam o iRobertsI (Il999h : Ijones and RobertsI (l2000l) : iBuffett 
(l20Q3h . 



2 Dynamo equtions 

2.1 Equations in physical space 

The geodynamo equations for the incompressible fluid (V • 
V = 0) in the volume of the scale L rotating with the angular 
velocity $7 in the Cartesian system of coordinates (x, y, z) in 
its traditional dimensionless form can be written as follows: 
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^ + (V.V)(T + To) = AT. 



The velocity V, magnetic field B, pressure P an d typica l 
diffusion time t are measured in units of k/L, y^2flKiJ,p;, 
pK^ /L^ and L^/k respectively, where k is thermal diffusiv- 

ity, p is density, /i permeability, Pr = — is the Prandtl num- 

V 

ber, E = ^ is the Ekman number, v is kinematic viscos- 

2VLI? 
ity, 77 is the magnetic diffusivity, and q = K/77 is the Roberts 

is the modified Rayleigh number, a 



number. Ra = 



2Q.K 



is the coefficient of volume e xpansion, bT is the unit of tem- 
perature, see for more details I JonesI (120001) . g^ is the gravita- 
tional acceleration, and Tq = 1 — 2; is the heating below. The 
problem is closed with the periodical boundary conditions in 
the (x, y) plane. In z-direction we use simplified conditions 



dCattaneo et al l2003h : Tq 
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2.2 Equations in wave space 

To solve the problem ([T]) we apply pseudo-spectral approach 
( Orszag. 1971) fre q uently used in geodynamo simulations 
Jjones and RobertsI I2OOOI: iBuffettL l2003h . Equations are 
solved in the wave space. To calculate the non-linear terms 
one needs make the inverse Fourier transform, then calculate 
a product in the physical space, make a Fourier transform of 
the product, and finally calculate a derivatives in the wave 
space. After elimination of the pressure using conditions of 



the free divergency k-V = 0, k-B==Owe come to 
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Vk = - 



k-Fk 
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k = kjskp^ /3 = 1 , 



Pr^^V X (V X V)+Rari, 



1^ X V + (B- V)B 



(3) 



For integration in time we use explicit Adams-Bashforth 
(AB2) scheme for the non-linear terms. The linear terms are 
treated using the Crank-Nicolson (CN) scheme. To resolve 
the diffusion terms we use the known trick which helps to 
increase the time step significantly. Consider equation 
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rewrite it in the form 



dt 
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(5) 



and then the CN scheme is applied. 

The most time consuming part of our MPl code is EFT 
transforms. To make our code more efficient we use vari- 
ous modifications of known EFT transforms which take into 
account special kinds of the symmetry of the fields. The op- 
timal number of processors for the grids 128^ is n ^ 50. The 
scalability tests demonstrates even presence of the superac- 
celeration for the number of processors < n. 



3 Basic properties of the fields 

Here we consider simulations without rotation similar to 



Meneguzzi and Pouqueti (119891) and with rotation for two 



regimes with different amplitudes of the heat sources: 

NR: Regime without rotation, Ra = 6 • lO*', Pr = 1, E = 1, 
q = 10, Re - 2.5 • 10^. 

Rl: Regime with rotation, Ra = 1.3 • 10^, Pr = 1, E = 
2- 10-*^, q= 10, Re- 1.6-103. 

R2: Regime with rotation, Ra = 2.1 • 10^, Pr == 1, E = 
2 . 10-6, q = 10, Re ~ 3 • 10^. 
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Fig. 1. Regime Rl. Sections of temperature distribution. . All sections correspond to the middle of the cube. The field ranges are (0, 1) 
left, (0.46, 0.54) - right. 
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Fig. 2. Distribution of \4-component of the velocity field with ranges (—248, 253), ( — 143, 144). 



The first (NR) regime is clos e to the typical Kolmogorov con - 
vection, see in more details iMeneguzzi and Pouqueti (Il989l) . 
Including of rotation (regime Rl) Fig. 1-5 leads to transform 
of the isotropic convective structures to th e cyclonic state 



with horizontal scale ^ E ' (kc ^ E ' ) IChandrasekhar 



l/3^ 



( 119611) . Including of the magnetic field (the full dynamo 
regime with magnetic energy comparable with the kinetic en- 
ergy on the order of magnitude) d oes no t chang e structure of 
the convective patterns too much ljonesi (|200Cr). In the same 



time spectra of the magnetic energy are quite different and 
have no well-pronounced maximum at kc. 



Increase of Rayleigh number leads to decrease of the rel- 
ative role of rotation and should decrease peak at the kinetic 
spectra energy what is accordance with spectra for regime 
R2, Fig. 5. In principle further increase of Ra should lead to 
the original Kolmogorov state, similar to NR with spectrum 
law ^ k~^/^ Fig. 5. However we emphasize, that informa- 
tion on the spectra is not enough to judge if the role of ro- 
tation is negligible or not and additional analysis is needed. 
The argument is the follows: rotation leads to degeneration of 
the third dimension (along z-axis) Batchelor (1953) . On the 
other hand in isotropic two-dimensional systems spectraum 
of the kinetic energy also have — 5/3-slope, however direc- 
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Fig. 3. Distribution of Vi-component of the velocity field with ranges (—675, 701), (—153, 157) 
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Fig. 4. DistributionofBz -component of the velocity field with ranges ( — 1.44, 1,14), (—1,88, 2.37) 



tion of the energy transfer in the system is inverse. In con- 
trast to the three dimensional turbulence where energy trans- 
fers from the small wave number, where energy is injected, 
to larger dissipative wave number, in two dimensional turbu- 
lence energy transfers from the large wave numbers to the 
small oneo So far the quasi-geostrophic turbulence inherits 
prop erties of the bo t h systems 2D and 3D, the inverse cas- 
cade |HossaiiJ(ll994|);|Constantin|(12002|), we plan to consider 
behavior of energy fluxes in the wave space more carefully. 



4 Energy fluxes 



To analyze e nergy transfer in the wave space we follow 
FrischI (Il995h , Decompose physical field / in sum of low- 
frequency and high-frequency counterparts: /(r) = /<(r) + 
/^(r), where 



There is also a direct cascade of enstrophy. 



/<(r)= E /^e 



ikr 



\k\<K 



/>(r)= E /fc^' 

\k\>K 



kr 



(6) 
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Fig. 5. On the left is spectra of the kinetic energy for NR (red), Rl (green), R2 (blue). On the right is a spectra of the magnetic energy. The 
straight line corresponds to the Kolmogorov's spectrum ~ fe~^". 
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Fig. 6. The normalized fluxes of kinetic Tk (on the left) and magnetic Tm (on the right) energies in the wave space: NR (red), Rl (green), 
R2 (blue). 



For any periodical / and g one has relation lFrischI ( 119951) : 
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<!{v)>=V'^ j f{v)dv^ 



(7) 
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means averaging of / over the volume V. Multiplying the 
Navier-Stokes equation by V< and induction equation by 
B< lead to equations of the integral fluxes of kinetic Ek = 
V'^ /2 and magnetic Em = B^/2 energies from k > K to 
k < K: 

Ha' =< (V X rotV) • V< >, 
Hm =< rot (V X B) • B< > 



(9) 



and for the flux of the Lorentz work: 

Ul =< (rotB X B) • V< > . 
Introducing 

dUK 



TK{k) = -- 



dk 



leads to obvious relation for Ek in fc-space: 
dEK{k) 



dt 



^T(k) + F{k)+D{k), 



(10) 



(11) 



(12) 



where k = |k|, T{k) is a flux of the energy from harmonics 
with different fc, F{k) is a work of the external forces and 
D{k) = —k'^Exik) is a dissipation. The accurate form of T 
is: 



Tk^ 
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Fig. 7. The normalized iluxes of the generation term ~Tl (on the left) and advective c Tn (on the right) terms in the wave space: NR (red), 
Rl (green), R2 (blue). 



Taking into account that rot (V x B) = 
(B ■ V) V one has: Tm ^ Tn - Tl, where 



(V-V)B 



Tn = — 



an 



N 



dk 



TN{k)dk = Q, 



Hat =< ((V-V)B<)-B< > 



fc=0 



(14) 



Fig. 6 presents fluxes of kinetic Tk and magnetic Tm en- 
ergies for the mentioned above regimes. Regime NR for 
Tk demonstrates well-known behaviour for the direct Kol- 
mogorov's cascade in 3D. For the large scales T^ < 0, 
these scales are donors, they provide energy to the sys- 
tem. On the other hand harmonics at the large k ab- 
sorb energy. The two-dimensional turbulence exhibits 
mirror-symmetrical behaviour relative to the axis of absciss 
iKraichnan and Montgomery ( 1980 ). In this case the energy 
cascade is inverse. 

Rotation changes behaviour of fluxes of kinetic energy es- 
sentially. The leading order wave number is fc^ For k > kc 
we also observe the direct cascade of energy Tk > 0. The 
maximum of Tk is shifted relative to the maximum of the 
energy to the large k as larger as larger Re. For k < k^ 
behaviour is more complex: for the small k the inverse cas- 
cade of the kinetic energy takes place Tk > 0. On the other 
hand for the larger region of fc (0 . . . fee) we still have the di- 
rect cascade Tk < 0. Increase of Re leads to narrowing of 
the region with inverse cascade and increase of inverse flux. 
One can suggest, that change of the sign of the flux Tk at 
k < kcis connected with appearance of the non-local energy 
transfer: so that energy to the large-scales ki comes from 
modes jkal - Ikgl > |ki|, ki = ka + kgiWaleffe (1992). 
In absence of the magnetic field maximum of TK{k = 1) 
appears. So, in case with rotation there are two cascades of 
kinetic energy (direct and inverse) tale place simultaneously. 

Now we consider the magnetic part. In contrast to Tk, Tm 
includes not only advective term but the generative term as 
well. This leads to positiveness of the integral Tm over all k. 



Moreover Tm is positive for any k. Position of maximum of 
Tm is close to that ones in spectra of Em , Tk ■ 

It is evident, that for the planetary cores distance between 
maxima in fluxes Tm for NR and Rl, R2 can be quite large, 
however not so large as kc- This statement concerned with 
condition on magnetic field generation which holds when the 
local magnetic Reynolds number Tm > 1 at the scale 1/fc: 
> 1 and that for the planets r/ ;» v. In the same 



kr] 



moment fluxes at the small k are small, i.e. system is in the 
state of the statistical equilibrium: dissipation at the small 
scales is negligible. 

Now we examine the origin of the magnetic energy at the 
scale 1/fc: does it concerned with the energy transfer from 
the other scales either it is a product of real generation at this 
scale? 

Fig. 7 demonstrates fluxes of — T^, concerned with mag- 
netic field generation. The maximum of generation term 
without rotation is at the large scale, while for the rotating 
system it is at ^ 1/fcc- Interestingly, that for the rotating sys- 
tem there is a region —Tl < for the large fc, where mag- 
netic field reinforce convection. For regime NR — T^ drops 
quickly because of the kinetic energy decrease Fig. 5. As a 
result we have: for the rotating system magnetic field is pro- 
duced by the cyclones, while for the non-rotating system the 
large-scale dynamo operates. 

Now we estimate the role of the advective term T/v sepa- 
rately. For the non-rotating Tjv and Tk are similar: the direct 
cascade takes place. For rotating system region fc ~ fee is a 
source of energy. In contrast to Tk, Tn has not positive re- 
gions at the small fc, i.e. the inverse cascade of the magnetic 
energy, concerned with the advective term in this region is 
absent. We draw attention to the amplitudes of the fluxes 



Tm, ~Tl, Tm'- for aU three cases hold 



T 



Ml 



Tr, 



10" 



i.e. 



two fluxes of the magnetic energy in the wave space with 
opposite directions exist. The first flux concerned with the 
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traditional energy transfer of the energy over the spectrum 
(advective term) and the flux of the Lorentz work. The re- 
gions of the maximal magnetic field generation coincide with 
the regions of the most effective magnetic energy transfer Tjv 
(from small k to the large k). Such a balance leads to equipar- 
tition state, when dissipation takes place on the large k. 

Note, that the full magnetic flux Fig. 6 Tm is localized at 
fc ^ 1 . For the non-rotating system it is because th e mean 



helicity and a-effect are zero Zeldovich et al.l (1198 3). Thus 



the inverse cascade of the magnetic energy to the small k is 
absent. 

For the rotating system here is a balance of the energy in- 
jection due to the Lorentz force and its wash-out because of 
advection. The latest effect reduces a-effect. 



5 Conclusions 

The magnetic fields of the planets are the main sources of 
information on the processes in the liquid cores at the times 
10^ — lO'^yy and more. If the poloidal part of the magnetic 
field is observable at the planet surface, then the largest com- 
ponent of field (toroidal) as well as the kinetic energy distru- 
bution over the scales is absolutely invisible for the observer 
out of the core. Moreover, due to finite conductivity of the 
mantle even the poloidal part of the magnetic field is cutted 
off at fc ^ fcc. In other words the observable part of the 
spectra at the planets surface is only a small part (not even 
the largest) of the whole spectra of the field. That is why 
importance of the numerical simulation is difficult to overes- 
timate. Here we showed that in considered quasigeostrophic 
state the both cascades (direct and inverse) exsist simulta- 
neously. That is a challenge for a turbulent models in the geo- 
dynamo. The other interesting point is a balance of the mag- 
netic energy flux due two mechanisms: advection and gen- 
eration. This balance is reason why the exponential growth 
of the magnetic energy is stopped at the end of the kinematic 
dynamo regime. 
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